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We apply the Jansson-Farrar JF12 magnetic field configuration in the context of point source 
searches by correlating the Telescope Array ultra-high energy cosmic ray data and the IceCube- 
40 neutrino candidates, as well as other magnetic field hypotheses. Our field hypotheses are: no 
magnetic field, the JF12 field considering only the regular component, the JF12 full magnetic field, 
which is a combination of regular and random field components, and the standard turbulent magnetic 
field used in previous correlation analyses. As expected from a neutrino sample such as IceCube-40, 
consistent with atmospheric neutrinos, we have found no significant correlation signal in all the cases. 
Therefore, this paper is mainly devoted to the comparison of the effect of the different magnetic 
field hypotheses on the minimum neutrino source flux strength required for a 5 ct discovery and the 
derived 90% CL upper limits. We also incorporate in our comparison the cases of different power 
law indices a = 2, a = 2.3 for the neutrino point source flux. Our results indicate that the discovery 
potential for a point source search is sensitive to changes in the magnetic field assumptions, being 
the difference between the models mentioned before between 15% and 26%. Finally, as a collateral 
result, we implement a novel parameterisation of the JF12 random component. 


I. INTRODUCTION 

One of the most important quests in particle astro¬ 
physics is to identify the sources that produce ultra-high 
energy cosmic rays (UHECR). One way to achieve this 
goal is to directionally correlate UHECR with astrophys- 
ical neutrinos. This correlation would not only indicate 
us the sites of hadronic acceleration, but also their non 
single-shot transient nature. On the other hand, a posi¬ 
tive correlation signal requires sources producing a simi¬ 
lar order of fluxes of UHECR (e.g.protons) and neutrinos. 
This condition is fulfilled by sources with a proton inter¬ 
action opacity r > 1. Meanwhile, sources with r 1 
or r ;g> 1 would either produce almost exclusively pro¬ 
tons with a small associated neutrino flux or produce 
only neutrinos, absorbing the corresponding protons, re¬ 
spectively. A mean free path length can be linked to the 
interaction between the protons escaping from the source 
and the Cosmic Microwave Background (CMB), defining 
the so called GZK sphere of radius ^ 100 Mpc [TJ [2]. 
This would give us the limit of the farthest distance of 
observable UHECR sources, an implicit condition for any 
correlation analysis. 

One key ingredient in a directional correlation analy¬ 
sis between UHECR and neutrinos is to have an accu¬ 
rate determination of the UHECR magnetic deflection. 
Since the UHECR are charged particles, in their travel 
to Earth, they are going to be deflected due to its in¬ 
teraction with the magnetic field inside and outside the 
Galaxy, the latter being true if they are coming from an 
extragalactic source. 

This type of analysis has been perfomed under different 
magnetic field hypotheses. Examples include correlation 
studies between UHECR from the Pierre Auger Obser¬ 
vatory (PAO) and neutrino events from the ANTARES 
Telescope [3] and between UHECR from Telescope Array 
(TA) and neutrino events from IceCube |1]. In all these 


analyses, an energy independent magnetic deflection has 
been used. This kind of estimation is reasonable since 
we are not certain about the description of the Galac¬ 
tic or the Extragalactic magnetic field. However, it is 
interesting to examine the impact on correlation analy¬ 
ses, calculating the magnetic deflections under different 
Galactic Magnetic field models, taking into account en¬ 
ergy dependent considerations in their estimations. Fol¬ 
lowing this idea, a study of cosmic ray deflections from 
Centaurus A was proposed in in which the more re¬ 
cent version of the Jansson and Farrar 0[7] (JF12) mag¬ 
netic field model has been implemented. The JF12 is an 
improved model of the Galactic magnetic field, includ¬ 
ing a regular and random field component. This model 
fits very well the WMAP7 Galactic synchrotron emission 
map and more than 40,000 extragalactic Faraday rota¬ 
tion measurements. The JF12 model also includes an 
out-of-plane component and striated-random fields. 

In this work we use four different magnetic field hy¬ 
potheses: no field, the regular component of the JF12, 
the full JF12 field model, and the standard turbulent 
magnetic field. These hypotheses are initially used to 
study the correlation between the Telescope Array ultra- 
high energy cosmic ray data [8] and the IceCube-40 neu¬ 
trino candidates [5]. Since no correlation is expected, the 
main analysis of this paper is the comparison of the neu¬ 
trino flux requirement for a 5tT discovery and the 90% CL 
upper limits, under the four magnetic field hypotheses. 
The paper is divided as follows: in section 2 we describe 
the magnetic field models used in this study and we also 
present the parametrization that we have developed of 
the angular deflections from the random component of 
the JF12 field. In sections 3 and 4 we outline the method 
used to find correlations between UHECR and neutrinos 
and we describe how to include magnetic field deflections 
in our statistical analysis. In section 5 we present our 
results and in section 6 our conclusions. 
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II. GALACTIC MAGNETIC FIELD MODELS 


The angular deflection S of an ultrarelativistic particle 
of charge Ze and energy E due to the magnetic field is 
described by: 


5 oc 


Ze 

~E 


p X Bds 


( 1 ) 


This formula reflects the inverse proportionality between 
S and E which is a key issue for understanding the pa- 
rameterisation presented in section 3. 



Angular deflection in degrees 



significant. In fact, this is particularly significant if the 
cosmic rays are composed by heavy nuclei such as iron 
and may be deflected by a few tens of degrees even for 
energies above 10^° eV. 

For distances up to 500 Mpc and energies E > ix 10^® 
eV, typical deflections from EGMF are smaller than the 
resolution of UHECR detectors (see, e.g., dH])- Thus, if 
the UHECR sources are within the GZK sphere of radius 
~ 100 Mpc, we may neglect EGMF deflection for energies 
E > 3 X 10^® eV, which is the case in our analysis. 

Our correlation analysis is mainly tested under two 
different Galactic magnetic held models. One of these 
models is the JF12 GMF model [6l [7], designed to ht 
the rotation measures and polarized synchrotron data, 
as we have already mentioned before. The JF12 model 
separates the magnetic held into three components: reg¬ 
ular (coherent), striated and purely random helds (ran¬ 
dom components). All components extend up to r = 20 
kpc, where r is the in-plane radius with origin at the 
galactic center. The regular component (or large-scale 
held) is a superposition of three helds: spiral disk held, 
toroidal halo held and the X-shaped poloidal halo held. 
Figurej^shows the dehection from the coherent held for a 
backtracked 57 EeV proton. Dehections are stronger for 
backtracked directions below the Galactic plane and are 
strongest when the particle propagates across the galac¬ 
tic center. 

In addition to the regular component of the JF12, the 
striated random held is fully aligned to all components 
of the regular held. The hnal component is an isotropic 
random small-scale held. The coherence length of these 
helds is expected to be of order 100 pc or less. The GMF 
will be modelled using the best-ht parameters in [i [7] 
and assuming a coherence length of 60 pc for the fully 
random component. 

The other model, which we denote by “standard”, was 
used in Ref. [3112] for similar UHECR-neutrino correla¬ 
tion studies. In this model the dehections follow a Gaus¬ 
sian distribution with width 3°. In |13] , this method was 
changed to a Fisher distribution with an energy depen¬ 
dent concentration parameter. 


FIG. 1. Top panel: Angular deflection from the regular com¬ 
ponent of the JF12 Field for a backtracked proton of energy 
E = 57 EeV. The galactic plane is indicated by the thick 
black line. Map drawn using HEALPix [14]. Bottom panel: 
Deflection of PAO and TA UHECR events under the regular 
component of the JF12 held. The shaded region marks the 
TA hotspot. Maps in equatorial coordinates. 

The dehection of a charged particle as it trav¬ 
els through a turbulent extragalactic magnetic held 
(EGMF), with strength of order 6nG, is proportional to 
'/D\Z\/E where D is the source distance, Z the charge 
number of the UHECR and E its energy [TU] [H] . These 
dehections are typically < 1° for protons with energy 
above 40 EeV and a propagation distance of order ^ 1 
Gpc. Meanwhile, dehections caused by the galactic mag¬ 
netic held (GMF), with strength of order 6^G, may be 


A. Random component parameterization of the 
JF12 GMF 

A secondary contribution of this work is the implemen¬ 
tation of a parameterization of the random component 
of the JF12 GMF model. This analytic expression would 
help us in having a better insight into the behaviour of 
this component while saving computational time in the 
calculations of our analysis. This parametrization is valid 
for small dehections, a condition that restricts the energy 
range of UHECR, which are protons, to energies greater 
than 57 EeV. We have chosen protons with this energy 
since the dehections are going to be small, which is a 
requirement for our results to be valid. In addition, for 
the backtracking method, we have used CRPropa 3 [TS] 
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FIG. 2. Left panel: R.M.S. deflection due to random field components for 100 EeV protons, using our parameterization (Eq. 
(5)). Right panel: The same plot by backtracking for several realizations of the full JF12 Field (simulation). Map drawn using 
HEALPix [II]. 
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FIG. 3. Deflection of the 69 UHEGR reported by PAO in 
m and the 72 TA events in Galactic coordinates. The blue 
crosses are the arrival directions of the cosmic rays as seen 
from Earth. The red dots correspond to the backtracked 
events using the regular component of the JF12 Field and 
each event is connected to its observed arrival direction by 
a black line. The gray circles have radii 2.0 hr.m.s (see Eq. 
(5)), which correspond to the additional deflection from the 
random components of the JE12 Field. 


to propagate cosmic rays under the JF12 magnetic field 
model, using the best fit values given in |^, for the co¬ 
herent field parameters, and in [7], for the random field 
parameters. 

Our procedure for obtaining the parameterisation goes 
as follows: firstly, for a given initial arrival direction (ob¬ 
served at the Earth) Pq = {lo,bo) of a proton with en¬ 
ergy E, where {I, b) are in the galactic coordinate sys¬ 
tem with —180° < I < 180° the galactic longitude and 
—90° < b < 90° the galactic latitude, we backtrack it 
through the coherent field to the position Pc = {lc,bc)- 
The latter is the final position reached by the proton 
before it leaves the GMF. Secondly, starting again from 


Pq, we backtrack another proton with the same energy 
E through the full magnetic field (coherent plus random) 
until it arrives to its final position P = {l,b), before it 
exits the GMF. Then the angular distance 6 between Pc 
and P is calculated. After repeating this process over dif¬ 
ferent realizations of the random component of the field, 
we may construct a distribution of the angular distance 
( 5 . 

To functionally describe the deflection, we have fol¬ 
lowed [13] where the random field deflections are given 
by a Fisher distribution. Under the small deflection angle 
hypothesis, the probability distribution function / of the 
random field deflections can be described by a Rayleigh 
distribution: 

/(d) = 2A(5exp(-A(52) (2) 

where A = X{E)\ig^i,c is a fit parameter for a given lo,bo. 
We find that the optimal energy dependence of A is of 
the form 

= ^l\lo,bo^ + ^2\lo,boE'^ (3) 

This energy dependence was tested for energies 57 EeV< 
E <200 EeV by taking 15 energies logarithmically spaced 
in this interval. The parameters Ai , A 2 are given on 
a 50 X 50 grid distributed uniformly across —180° < 
lo < 180° and —0.9999 < sin(&o) < 0.9999. Using the 
2500 tabulated points as a sample, the weighted aver- 


ages (weight Wi = 

cos(6i)) are given by 


(Ai)^ 

^ (0.1 ±3.0) X lO-^EeV^ 



^ (2.6 ± 1.9) X 10"^EeV"^ 

(4) 

The distribution / also satisfies 


(d)= ^ 

dr.m.s. = ^/W) = ^ 

(5) 
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From Eq. (5) we see that the ^r.m.s deflection decreases 
for higher energies being that for energies above 100 EeV 
it goes like A « A 2 E^. The behaviour ^r.m.s fx E~^ is 
therefore recovered in the higher energy regime, in agree¬ 
ment with Eq. (1). In fact, this is the reason why there 
are no terms of order higher than E'^ in the expression 
for A. In order to test our proposed parameterization, we 
display in Eigurej^a comparison between this one and the 
numerical simulations, which are based on Runge-Kutta 
methods for the dr.m.s. of 100 EeV protons. These simula¬ 
tions are performed over 640 realizations of the full JE12 
field and with the sky divided into only 49152 arrival di¬ 
rections. It is clear the excellent agreement between our 
parameterization and the numerical calculations. 

Using / we can reduce the computational time by up to 
two orders of magnitude, in comparison with the numer¬ 
ical calculations. The function / only gives information 
on the angular deflection S, leaving for each 6 an allowed 
region, an annulus, for the final backtracked position. We 
assume a symmetrical distribution of the particles in the 
annulus. 

In Figure we show the effect of Galactic magnetic 
deflections on the 69 PAO UHECR events and 72 TA 
events. These deflections have been obtained due to the 
regular component, that produce a fixed separation, and 
the random component, depicted using our parameteri¬ 
zation. It is evident how the greatest deflection given by 
the random component occurs mainly in the vicinity of 
the galactic plane. 


III. GENERAL SCHEME OF THE 
CORRELATION ANALYSIS 

In order to calculate the correlations between IceCube 
neutrinos and TA UHECRs, we have followed the source 
stacking method outlined in [3- This method adds up 
the flux intensities from a group of single sources con¬ 
centrated in a small region of the sky. For the UHECR, 
we consider 72 TA UHECR events with energies E > 57 
EeV using Surface Detector data collected between May 
2008 and May 2013 [8]. We have chosen the energy cut 
at 57 EeV, since the backtracking method, in which our 
analysis relies on, is not able to deal with the large GMF 
deflections produced for energies lower than this limit. 
For the neutrino sample, we have used 12877 upward 
going muon neutrino candidate events recorded by Ice- 
Cube in the 40-string configuration between April 2008 
and May 2009 with a livetime of 375.5 days [9]. While 
we are aware that there are more recent neutrino sam¬ 
ples (e.g. m), the arrival direction of individual events 
are not available, providing a distribution in zenith angle 
bins instead. 

For the purposes of this study we also work in Equa¬ 
torial coordinates {a,S), where 0° < a < 360° is the 
right ascension and —90° < S < 90° is the declination. 
Given the geographical location of the IceCube neutrino 
observatory, the declination can be easily related to the 


zenith angle 9 via the relation i5 = 90 — 0. For the galac¬ 
tic magnetic field scenario, we study four cases: no field, 
only with the regular JF12 field, the full JF12 held (regu¬ 
lar plus random components) and the standard turbulent 
held. 

In order to give a complete picture of our correlation 
analysis it is important to comment on the differences 
between the study in [3] and ours. First, in contrast with 
the magnetic dehections used in (2, which are isotropic 
and independent of CR energy E, we are obtaining de¬ 
hections as a function of E and the arrival direction. The 
latter is an important difference since the dehections ob¬ 
tained via the backtracking method depend on the mag¬ 
netic held traversed by the cosmic ray. Second, we are in¬ 
cluding the dehections from the regular held components 
that we use for shifting the cosmic ray arrival positions, 
prior to using them for our correlation analysis. Third, 
we are considering in our study a variety of GMF models. 

A. SIGNAL AND BACKGROUND 

The expected number of signal neutrinos jls is given 
by: 

jls=T j dVt J dE^AeffiE^,cos6)^i,{E^) (6) 

where T is the livetime, A^ff is the IceCube muon- 
neutrino ehective area in the 40-string conhguration, E^ 
is the neutrino energy, 6 is the zenith angle, $ is the 
muon-neutrino source hux and dO, is the differential solid 
angle element. Effective areas are averaged over 30° 
zenith bins as presented in [9]. Our hypotheses for the 
neutrino hux are (X E~‘^ and (X based on 

recent IceCube results [TB]. For a better understanding 
of our results it is important to remark that the /is is 
implicitly dependent on declinations since they are de¬ 
pendent on the zenith angle. 

To obtain the background we generate 10® pseudo¬ 
experiments by scrambling the right ascension in the neu¬ 
trino data. Centered at each cosmic ray, we dehne a re¬ 
gion of points within an angular distance di that we call 
angular bins of size 4*. As mentioned before, we include 
dehections from the regular held components into the cos¬ 
mic ray positions by shifting them accordingly. For each 
the expected count /rb('I') of neutrinos and its stan¬ 
dard deviation ab are determined assuming a Gaussian 
distribution. 

The value of 4' yielding the lowest 90% Feldman- 
Cousins conhdence level (CL) mean upper limit (/igo) 
m is calculated, assuming background only. The value 
of (/Tgo) depends on /i(, only. We then calculate the ex¬ 
pected signal jis based on Eq. (6), but adding the ef¬ 
fects of angular resolution of the experiments by spread¬ 
ing the neutrino source according to a Gaussian distri¬ 
bution with standard deviation 2°. We also include an 
additional spread when dealing with the standard held 
and the JFI2 random held components. We denote by 
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/is the expected signal which embodies both effects (an¬ 
gular resolution and random magnetic field components). 
Because our magnetic field models are fixed before apply¬ 
ing the analysis, there are no associated trial factors in 
the calculations. 

We find the angular bin size 'h mrf which minimises 
the model rejection factor MRF=(/igo)//is- For this 
'^MRF, we determine the 5a 90% CL discovery poten¬ 
tial /is°’, calculated from the requirement P{n > fib + 
5ab\fi = fib + fi^'^) = 0.5, where n is number of events. 
We then find the normalisation constant C = 
which gives the strength of the source, from the equa¬ 
tion /ts*^ = fia{C) and solving for C. Note that while 
fF’f is the expected neutrino count of signal neutrinos 
summed over all sources, the value of C is the strength 
per source. 


malisation constant is similar in both the B = 0 

and reg. JF12 cases because of the negligible change 
in mrf - When we compare the full JF12 against the 
standard scenario, we see that /t®'^ is stronger in the full 
JF12 but is weaker. This odd behaviour comes 

from the different treatment of the magnetic deflections 
in these models. As mentioned before, the standard 
scenario assigns an energy independent deflection, with 
<5r.m.s. = 3.0°, which tends to overestimate the random 
deflections in comparison with the full JF12 (see Figure 
[^. As a consequence, for the standard case, fewer neu¬ 
trino events are enclosed within an angular bin size 'k 
compared to the full JF12, thus the normalisation con¬ 
stant needs to be higher in order to reach a given 

value of fi^/. A similar relation is seen in Table 2. 


IV. RESULTS 

We find that for the four GMF model scenarios: no 
field (B=0), only with the regular JF12 field (reg. JF12), 
the full JF12 field consisting of the regular plus random 
components (full JF12) and the standard turbulent field, 
there is no 5tT discovery after unblinding the data, be¬ 
ing all the observed neutrino events consistent with the 
background. This result is expected since the chosen neu¬ 
trino sample is consistent with atmospheric neutrinos and 
the compatibility between data and background is clearly 
shown in Figure]^ 

In this way, introducing the full JF12 field does not 
improve the intensity of a correlation signal between 
UHECR and neutrinos in comparison to other field as¬ 
sumptions. However, and in spite of this result, it is 
interesting to foresee the impact of choosing a particu¬ 
lar GMF model in the values required to satisfy a pos¬ 
itive correlation analysis. For that reason we consider 
worthwhile to extend our study in order to include the 
values of fi^f and the 90% GL upper limits of the nor¬ 
malisation constant, which we obtain from the relation 
M 90 (Mbj''^obs) = Ms(C); where /rgg is calculated according 
to the Feldman-Gousins prescription. 

In Tables I and 2 we show, for each of the GMF hy¬ 
potheses, the values of '^mrf, and the correspond¬ 
ing flux normalization constants E“<I>®°’ and of 

/ig*^ and 90% CL upper limits respectively. Tables 1 
and 2 were obtained by using the 72 TA events and 
$ 1 . oc E“^/<I>y oc E~^-^ signal respectively. 

According to each magnetic field scenario, the be¬ 
haviour of mrf and /t®°’ are similar in the sense that, as 
'i'MRF increases, so does fi^'^. This can be explained by 
the fact that fib and Uf, increase with increasing 'i’MRF, 
requiring a larger /i®°" to get the desired 5cr excess above 
the expected background. The reason why we get a 
slightly higher /tg°' for E = 0 than the reg. JF12 case 
is that the change in declinations, after performing the 
magnetic field deflections for the reg. JF12 component, 
turns out in a variation of A^ff and thus in /tg°'. The nor- 


500 



B=0 regJF12 Full JF12 Standard 
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FIG. 4. Distribution of expected background, observed events 
and discovery potential for different magnetic field hypothesis. 

It is also relevant to mention the quantitative differ¬ 
ences between the full JF12 and standard field models. 
In Tables 1 and 2, is smaller in the full JF12 by 

15% and 17% respectively, indicating that the strength of 
the point source required for a 5a discovery is somewhat 
overestimated in the case of the standard field. 

Regarding the dependence of the power law index a 
on our results, we observe that the relative differences 
between the normalization constants for different field 
models are independent of a. For the normalization con¬ 
stant itself we get higher values when a = 2.3. This can 
be explained by the steeply falling flux in this scheme, su- 
pressing contributions from the highest energy neutrinos 
at the tail, which have the largest effective area and are 
the most relevant for fig. The values of fi^'^ remain almost 
intact after a change in a, since the effect of this change 
is subdominant in our procedure. The results presented 
in Tables 1 and 2 are also shown in Figure which is 
valid because of the negligible differences in 'Lmaf be¬ 
tween both tables. The differences in the upper limits is 
justified by the deviations of the observed numbered of 
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events with respect to its expected background. If we ob¬ 
serve a large excess above the background, a larger flux 
per source consistent with the observed data are allowed. 
For this reason, in the full JF12 scenario, we have an 
upper limit ~ 30% higher than the other field models. 

Finally, in Table 3 we explore the effects on due to 
variations of the JF12 field model parameters and com¬ 
pare them to the results of our parameterization, assum¬ 
ing an signal flux and using the 72 TA event data 
sample. We chose the parameters Bq and (3, where Bq 
represents the magnetic field strength of the random halo 
component and j3 modulates the strength of the random 
striated component with respect to the regular compo¬ 
nent via the relation B^^^^ = PB^^g. There is an excellent 
agreement, of about 1%, between our parameterization 
and the JF12 numerical simulations, taking its best fit 
parameters. We achieve differences of at most 5% after 
varying both parameters by one standard deviation. 

V. SUMMARY AND CONCLUSIONS 

We have made a correlation analysis between the TA 
UHECR events and the IceCube-40 muon neutrino candi¬ 
dates under different GMF model assumptions. In par¬ 
ticular, we have introduced the JFI2 model, which ex¬ 
presses our best knowledge of the GMF. We have not 
found any correlations between the samples, regardless of 
the magnetic field hypothesis and the signal flux power 
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TABLE I. 5(t 90% CL discovery potential for the 72 TA events {E^ ^ signal flux). and are given in 

GeV cm“^s“^sr“^. 


Field Model ^ 

MRF (deg) 


per source 

E^^l° per source 

B = 0 

1.13 

60.5 

5.05 X 10"® 

1.02 X 10"® 

Reg. JF12 

1.10 

59.7 

5.04 X 10"® 

1.54 X 10"® 

Full JF12 

1.79 

94.9 

6.17 X 10"® 

2.49 X 10"® 

Standard 

1.19 

63.6 

7.27 X 10"® 

1.93 X 10"® 


TABLE 11. Same as Table 1 for an E^, signal flux 


Field Model 'I' 

MRF (deg) 


per source 

per source 

B = 0 

1.14 

61.2 

6.85 X 10-® 

1.44 X 10"® 

Reg. JF12 

1.10 

59.7 

6.89 X 10"® 

2.10 X 10"® 

Full JF12 

1.81 

95.0 

8.34 X 10"® 

3.65 X 10"® 

Standard 

1.19 

63.6 

1.00 X 10"'^ 

2.64 X 10"® 


TABLE III. 5a 90% CL discovery potential for the 72 TA events {E^ ^ signal flux) under the Full JF12 Magnetic Field with 
varying field parameters 


Field Model 

'^MRF (deg) 


per source 

Our work 

1.79 

94.9 

6.17 X 10"® 

Best fit 

1.81 

96.6 

6.22 X 10"® 

P + a 

1.80 

95.3 

6.22 X 10"® 

P -a 

1.70 

91.0 

6.21 X 10"® 

Bo + a 

1.83 

96.9 

6.36 X 10"® 

Bo — a 

1.90 

100.1 

6.06 X 10"® 

P + a, Bo + a 

1.99 

104.7 

6.37 X 10"® 

b 

1 

o 

b 

1 

1.70 

90.9 

6.08 X 10"® 



